N8 8-18956 

A Superelement Component Dynamic Synthesis Method 


M. L. Soni 


A method is presented for coupling a broad class of component 
dynamic models in the manner of direct stiffness assembly and 
is implemented in a general matrix manipulation program. 


INTRODUCTION 

A number of methods have been developed in the past to accomplish the basic 
objectives of the method. These methods are surveyed in detail in [1]. Each of 
these methods are characterized by the way the component dynamics is input and 
coupled to adjacent components. Hurty [2,3] proposed that the connect degrees of 
freedom (DOF) of a component were fixed or had a zero displacement. He then 
partitioned the modes of the structure into rigid body modes, constraint modes, and 
normal modes. Craig and Bampton [4] proposed a simplification of Hurty' s fixed 
interface method by dividing component modes into only two groups: constraint modes 
and normal modes. Bamford [5] added attachment modes to improve the convergence of 
the method. 

Goldman [6] developed the free interface method, employing only rigid body 
modes and free-free normal modes in component dynamic representation. This 
technique eliminates the computation of static constraint modes, but their advantage 
is negated by the poor accuracy of the method. Hou [7] presented a variation of 
Goldman's free- interface method in which no distinction is made between rigid body 
modes and free-free normal modes. Gladwell [8] developed "branch mode analysis" by 
combining free interface and fixed interface methods to reduce the order of the 
stiffness and mass matrices for individual substructures. The reduction procedure 
requires the knowledge of topological arrangement and dynamics of all the components 
in the model. In order to account for adjacent components Benfield and Hruda [9] 
introduced inertia and stiffness loading of component interfaces. The use of loaded 
interface modes is shown to have superior convergence characterisics. 

MacNeal [10] introduced the use of hybrid modes and inertia relief modes for 
component mode synthesis. Hybrid modes are substructure normal modes computed with 
a combination of fixed and free boundary conditions. Inertia relief attachment 
modes are attachment modes for components with rigid body freedoms. MacNeal also 
included residual inertia and flexibility to approximate the static contribution of 
the truncated higher order modes of a component. Rubin [11] extended the residual 
flexibility approach for free interface method by introducing higher order 
corrections to account for the truncated modes. Klosterman [12] more fully 
developed the combined experimental and analytical method introduced by MacNeal. 
Hintz [13] discussed the implications of truncating various mode sets and developed 
guidelines for retaining accuracy with a reduced size model. 
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Recent research has centered on the ccmparision of the various methods. Baker 
[ 1 -4 ] , for example, compares the constrained and free-free approach using 
experimental techniques and also investigates using mass additive techniques and 
measured rotational DOF [15]. This investigation was motivated by a need to find 
the best method for rigidly connected flexible structures. In this connection, the 
constrained method produced the best result. Klosterman [1 6] has shown the free- 
free method to be accurate for relatively stiff structures connected with flexible 
elements. This supports Rubin's conclusion [11] that the free-free method is at 
least as accurate when residual effects are accounted for. These conclusions are 
intuitive because the type of boundary condition imposed in the analysis that best 
represents the boundary of the assembled structure provides the best accuracy in the 
modal synthesis. 

Meirovitch and Hale [16] have developed a generalized synthesis procedure by 
broadening the definition of the admissible functions proposed by Hurty [1]. This 
technique is applicable to both continuous and discrete structural models. The 
geometric compatibility conditions at connection interfaces are approximately 
enforced by the method of weighted residuals. 

The method due to Klosterman [12] has been implemented in an interactive 
computer code SYSTAN [17] and that due to Herting [18] is available in NASTRAN. The 
latter is the most general of the modal synthesis techniques. It allows retention 
of an arbitrary set of component normal modes, inertia relief modes, and all 
geometric coordinates at connection boundaries. Both the fixed- interface method of 
Craig and Bampton, and the MacNeal's residual flexibility method, are special cases 
of the general method. Other analyses presented in the literature based on modal 
synthesis techniques are not incorporated into general structural analysis codes. 

It is desirable to have a synthesis method to couple different types of 
component dynamic models in a setting such as that of the finite element method so 
as to be able to realize the best advantages of both the component synthesis methods 
and the convenience and generality of the finite element method. The need for a 
such a capability arises when different types of dynamic models are used to 
represent the various components of a structure. This variance is dictated by the 
need to improve accuracy of component dynamic representation, availability of a 
certain kind of data, etc. Herting's work [18] meets some of these objectives. The 
procedure presented herein permits a broader class of component models in the manner 
of direct stiffness finite element assembly and can be implemented in a general 
matrix manipulation program. 

A principal feature of the work developed here is the component dynamic model 
reduction procedure that leads to an exact and numerically stable synthesis. In 
order to affect component coupling, neither the specification of external coupling 
springs nor an user-specified selection of independent coordinate is required. 
Existing synthesis procedures suffer frcm these drawbacks. Component dynamic models 
considered include free-free normal modes with or without interface loading, up to 
second order stiffness and inertia connections accounting for the effect of modal 
truncation, fixed interface modes, and also the physical coordinate components. The 
model reduction procedure involves interior boundary coordinate transformations 
which explicitly retain connection interface displacement coordinates in the reduced 
component dynamic representations. Interior coordinates may include physical, 
modal, or any admissible coordinates. Components in this reduced form are termed 
"superel ament" because they are a generalization of the conventional finite elements 
of structural mechanics. The problem of component dynamic synthesis is then reduced 
to the assembly of the superel ament. The direct stiffness approach and all 
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subsequent processing operations of the finite element method are then applicable. 
The formulation, implementation and verification aspects of the superelement 
component dynamic synthesis method are presented in the following sections. 

SUPERELEMENT MODEL REDUCTION 

The procedures for reducing component dynamic models to a form involving 
physical DOFs of the connection interface nodes plus seme additional DOFs related to 
the interior nodes are derived in the following. The component dynamic models in 
this reduced form possess a structure similar to that of the displacement finite 
element method and are therefore termed superelements. Four types of component 
dynamic models are considered: (1) finite element model; (2) fixed interface modal 
model with static constraint modes; (3) free interface modal model with residual 
flexibility attachment modes; and (4) free interface model with normal modes. The 
last type may also include any general admissible shape vectors and corresponding 
dynamic matrices as long as certain requirements for matrix partitioning and 
invertibility are satisfied. A system may involve any combination of the above 
types of component dynamics model since the models are reduced to a common form 
before assembly. 

The finite element and the constrained interface modal models naturally 
contains the connection interface DOFS and are already in the required superelement 
form. Constrained interface modal model based superelement synthesis is treated in 
[19]. The free interface modal model with residual flexibility attachment modes is 
reduced to the superelement as follows. Expressing component displacements X in 
partitioned form can be witten as 
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where <{> and G are respectively the normal mode and residual flexibility attachment 

mode matrices, q^ and _q^ are the associated generalized coordinates, and the 

subscripts B and 0 refer to the boundary and other DOFs. To obtain the superelement 
reduction transformation, following Martinez [20] solve the lower partition of Eq. 
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substituting the above in Eq. (1) leads to the superelement form associated with the 
residual flexibility attachment model: 
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The generalized coordinates q^ are now the participation factors of the modified 
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in effect free interface normal mode partitions minus a set of constraint modes 
internal partitions and therefore represent modes constrained at boundary degrees of 
freedom. The superelement equations of motion are obtained by transforming the 
component modal equation of motion as 
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where Jg g = G gg H gg Ggg ,and m^ and are generalized mass and stiffness 

matrices respectively, Ggg and Hgg are the generalized stiffness and mass 

associated with the residual flexibility attachment modes. The modification of the 
component generalized mass and stiffness matrices by the contributions from the 
residual flexibility of the deleted modes is clearly seen in the above equation. 

In the event only a truncated set of free interface normal modes is used to 
represent component displacements, the superelement reduction is obtained from a 
partial inversion of the modal matrix as described in the following. Consider the 
free interface normal mode transformation 
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and (n B xm K )size partitions of the component modal 
matrix corresponding to the interior and boundary degrees of freedom Xq and Xg, 
respectively. Consider the lower portion of Eq. (5) 


where | KQ and ^ are (n Q xm K ) 


(5) 




( 6 ) 


and partition ^^g into a nonsingular invertible square matrix <f> KB1 
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which requires that the number of kept modes be greater than the number of boundary 
degrees of freedom. Solving the above equation for q^ gives 
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the generalized coordinates are called reduced modal coordinates since they are 

associated with modified modes as seen in the above equation. Combining transform 
Eqs. (6) and (9), a net transformation from physical coordinates X to free interface 
normal mode superelement coordinates can be written as 
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where and are Pactiticcs of <j>^Q created using the partitioning information 


generated in Eq. (8). Analogous to the Ritz transformation for constrained 
interface modal models, and the transformation of Eq. (4) for residual flexibility 
model modal, the left and right partitions of the transformation of Eq. (10) may be 
interpreted as expressing modal matrices of modified normal modes and boundary 
constraint modes, respectively. The component equations of motion in superelement 
coordinates (c^, Xg) can be obtained using the transform of Eq. (11) as 
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using physical coordinate component model, or 
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from modal coordinate model. In the above the generalized mass and stiffness 
matrices are used to include the case where the columns of the modal matrix in Eq. 
(12) are not orthogonal. (Mass and stiffness loaded component modes are 
nonorthogonal with respect to the original unloaded mass and stiffness matrices, for 
example.). Further, it is to be emphasized again that the reduced modal 
representation given in Eqs. (9) through (12) is obtainable from any given free 
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interface ’’modal" representation provided the partitioning conditions leading to 
Eq. (8) can be met. The columns of modal matrix in Eq. (5) need not be normal 
modes. The only necessary condition is that they be linearly independent, adequate 
in number, and be from a complete set, i.e., a linear combination of them should be 
capable of representing the deformation shapes of the component undergoing motion 
within the compound of the built-up system. 

SUPERELEMENT SYNTHESIS 


The reduction procedures described in the preceding section yield component 
models with a common characteristic: the boundary degrees of freedom are explicitly 
retained and the internal degrees of freedom are transformed to reduced generalized 
coordinates. The equations of motion of a superelement component can be expressed 
as 
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Each component substructure of a given built-up system can be expressed in the above 
format. Analogous to the matrix equation of a displacement based finite element 
model, the generalized coordinates may be interpreted as internal DOFs. Assuming 
that the global structure as partitioned into N components , the assembled global 
matrices M, C, K relating the global displacement vector X and the global force 
vector £ can now be synthesized by the direct stiffness approach fran the component 
matrices of Eq. (13). The direct stiffness approach postulates equal nodal 
displacement and nodal equilibrium of forces at the connection interfaces of 
adjacent components, thus 
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No compatibility conditions exist among the reduced generalized coordinates £ and 
remain intact. The coupled system coordinate vector thus becomes 
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so that the global mass, damping, and stiffness matrices are given by 
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where B is the portion of the component to global transformation matrix § 
corresponding to the a n component. 
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Y being the vector of uncoupled superelement coordinates. 
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and representing the compatibility conditions for the degrees of freedom at the 
boundaries of the component a. 


The built-up system mass, damping, and stiffness matrices are of the following 

form. 
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where the coefficients a) . are the mass, damping, or stiffness coefficients 
matrices given in Eq. (19), and A is the corresponding global matrix. The diagonal 
submatrices a'‘ ot ' in Eq. (19) correspond to the reduced generalized coordinates of 
the component^and are uncoupled frcm other components. 
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The submatrices are symmetric and in general fully populated. The system governing 
equations can be expressed as 


M X + C X + K X = f 


( 20 ) 


which itself is in a superelement form since the system basis coordinate X contains 
the boundary degrees of freedom explicitly. This is particularly useful frcm the 
standpoint of coupling the system of Eq. (20) to a higher level superelement. The 
recovery of component displacements frcm the system displacement vector X is 
simply a back substitution and transformation process. 


NUMERICAL EXAMPLE AND CONCLUSIONS 


The above formulation is implmented as a pre- and post-processor to a general 
matrix manipulation computer program. The following presents a sample synthesis 
problem involving free interface and residual flexibility modal models solved using 
the developed software to demonstrate its working. In Figure 1 the component A 
represents a lumped parameter model of an aircraft, while that labeled B represents 
a model of a store to be attached at the tip of the wing. The types of dynamic 
models employed for each component are as follows. For Aircraft: (1) free interface 
normal modes followed by a transformation to superelement coordinates, and (2) Free 
interface normal modes plus a residual flexibility attachment mode followed by a 
transformation to superelement coordinates. And for the Store: (1) free interface 
normal modes followed by superelement coordinate transformation, and (2) physical 
coordinate model. For comparison purposes the problem is also solved using existing 
MacNeal Method and Rubin Method of synthesis. The results of the synthesis are 
shown in Table 1 . The exact results are obtained by solving the eigenproblem of the 
built-up system without partitioning. The superelement synthesis using complete 
mode sets of the components leads to exact system synthesis as expected. 

Superelement method using severely truncated component mode set along with residual 
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flexibility attachment modes predicts system frequency with an accuracy better than 
that of the MacNeal method and equal to that of the Rubin's second order method. 

It must be noted that in this study different types of component dynamics 
models are synthesized by virtue of the superelement formulation, and that the 
accuracy of system synthesis is entirely governed by the type of component dynamic 
models and not by the coupling procedure. 


TABLE 1 DIRECT AND SYNTHESIZED SYSTEM NATURAL FREQUENCIES (RADS. /SEC) 
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EXACT 
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METHOD 

METHOD 



A 

B 

(B) 

(B) 

1 

0 . 

0 . 

0 . 

0 . 

0 . 

2 

52.417 

52.417 

52.418 

52.419 

52.418 

3 

69.348 

69.348 

69.829 

69.886 

69.829 

4 

72.818 

72.818 

291 .697 

504.67 

291 .697 

5 
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Fig. 1 Lumped Parameter Aircraft-Store System 
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